o 



^ 



> 






A NUMERICAL STUDY OF INFINITELY RENORMALIZABLE AREA-PRESERVING MAPS 

DENIS GAIDASHEV AND TOMAS JOHNSON 



Abstract. It has been shown in (Gaidashev and Johnson 2009b) and (Gaidashev et al 2011) that infinitely renormal- 
izable area-preserving maps admit invariant Cantor sets with a maximal Lyapunov exponent equal to zero. Further- 
more, the dynamics on these Cantor sets for any two infinitely renormalizable maps is conjugated by a transformation 
that extends to a differentiable function whose derivative is Holder continuous of exponent a > 0. 

In this paper we investigate numerically the specific value of a. We also present numerical evidence that the 
00 , normalized derivative cocycle with the base dynamics in the Cantor set is ergodic. Finally, we compute renormalization 

eigenvalues to a high accuracy to support a conjecture that the renormalization spectrum is real. 



1. Introduction 



00 

Q 

jrt \ Numerical investigations of universality properties in dynamics have been historically important, and, in general, 

motivated the renormalization approach to universality. Such studies were behind the pioneering discovery of the 
Feigenbaum-CouUet-Tresser period doubling universality in unimodal maps (Feigenbaum 1978), (Feigenbaum 1979), 
(Tresser and Coullet 1978). 

Universality problems are typically approached via renormalization. In a renormalization setting one introduces 
a renormalization operator on a functional space, and demonstrates that this operator has a hyperbolic fixed point. 
fvS ' Period-doubling renormalization for dissipative two-dimensional maps has been extensively studied in (Collet et 
■r^ . al 1980, de Carvalho et al 2005, Lyubich and Martens 2011). Compared to the Feigenbaum-CoUet- Tresser one- 
CO ' dimensional renormalization, the new striking feature of the two dimensional renormalization for highly dissipative 
r — maps is that the restriction of the dynamics to their Cantor attractors is not rigid. Indeed, if the average Jacobians 

C^ of F and G are different, for example, hp < bo, then the conjugacy F\cp ^ G\ca is not smooth, rather it is a Holder 

I I continuous function with a definite upper bound on the Holder exponent: a < 5(1 + ^ttF') *^ ^- ^^ turns out that 
IIV ' the period-doubling renormalization for area-preserving maps is very different from the dissipative case. 

A universal period-doubling cascade in families of area-preserving maps was observed by several authors in the 
early 80's (Derrida and Pomeau 1980, Helleman 1980, Benettin et al 1980, Bountis 1981, Collet et al 1981, Eckmann 
et al 1982). The existence of a hyperbolic fixed point for the period-doubling renormalization operator 



Rekw[F]=^Ap^oFoFoAf, 

where Ap{x,u) = (Xpx, fipu) is an _F-dependent linear change of coordinates, has been proved with computer- 
assistance in (Eckmann et al 1984), and in (Gaidashev et al 2011). The dynamics of the area- preserving fixed point 
map has been studied in (Gaidashev and Johnson 2009a, Gaidashev and Johnson 2009b, Johnson 2011, Gaidashev 
et 0^2011) 

We have proved in (Gaidashev and Johnson 2009b) that infinitely renormalizable maps in a neighborhood of the 
fixed point of (Eckmann et al 1984) admit a "stable" Cantor set, that is the set on which the maximal Lyapunov 
exponent is zero. We have also shown in the same publication that the conjugacy of stable dynamics is at least 
bi-Lipschitz on a submanifold of locally infinitely renormalizable maps of a finite codimension. This has been further 
improved in (Gaidashev et al 2011), where it is shown that the conjugacy of dynamics on the Cantor set in fact 
extends to a differentiable map whose derivative is Holder continuous with an exponent a > 0. This seems to be in 
stark contrast with the case of dissipative maps. 

In this paper we continue the study of the dynamics on the stable Cantor set, and we provide numerical evidence 
in support of the following two conjectures that will be described in detail in the paper: 
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i) the normalized derivative cocyle on the Cantor set is ergodic with respect to an absolutely continuous invariant 
measure on the circle; 

ii) the Holder exponent a satisfies a > 0.070347 for all infinitely renormalizable maps in a neighborhood of the fixed 
point. 

We also provide numerical evidence for the following conjecture 

Hi) the area-preserving period doubling operator has a real spectrum. 

The last conjecture can be compared with the corresponding conjecture in the one-dimensional case (Christiansen 
et al 1990). In particular, we compute approximations to the leading 100 eigenvalues of the operator, and they 
appear to be real, up to numerical accuracy. 

2. RENORMALIZATION for AREA-PRESERVING REVERSIBLE TWIST MAPS 

An "area-preserving map" will mean an exact symplectic diffeomorphism of a subset of M^ onto its image. 
A reversible area-preserving map that satisfies a twist condition 

duin,F{x,u))^0 

everywhere in its domain of definition can be uniquely specified by a generating function s: 

(1) f M-f ' 

\-s{y,x)J \s{x,y) 

The universality phenomenon in area-preserving maps can be explained rigorously if one shows that the following 
renormalization operator 

(2) Rekw[F]^K-/oFoFoKf, 

where A^ is some _F-dependent coordinate transformation, has a fixed point, and the derivative of this operator is 
hyperbolic at this fixed point. According to (Collet et al 1981) Kp can be chosen to be a linear diagonal transfor- 
mation: 

Kf{x, u) = {\x, ^u). 
The induced operator on generating functions will be denoted by TIekw- 

We will now summarize the definitions and theorems for the renormalization operator acting on generating func- 
tions originally due to the authors of (Eckmann et al 1982) and (Eckmann et al 1984). Part {i) — (iii) of Theorem 
I2.3l were proved in (Eckmann et al 1984), Theorem 12. 81 was proved in (Gaidashev and Johnson 2009b), and Theorem 
12.31 part {iv) and Theorem 12.91 were proved in (Gaidashev et al 2011). 

Consider the dyadic group, 

(3) {0,1}°° =^{0,1}", 

where lim stands for the inverse limit. An element w of the dyadic group can be represented as a formal power series 
w — >■ X]fc°=o ^fc+i^*^- The odometer^ or the adding machine, p : {0, 1}°° — > {0, 1}°° is the operation of adding 1 in this 
group. 

We will require a couple of more definitions for a precise statement of our results from (Gaidashev et al 2011). 

Definition 2.1. The Banach space of functions s{x, y) = ^^ -^q Cij{x — PYiy — /?)-', analytic on a bi-disk 

VpiP) = {(x, y)eC^:\x~P\< p, \y ^ I3\ < p}, 
for which the norm 

CXJ 

is finite, will be referred to as A^ {p) . 
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A^{p) will denote its symmetric subspace {s G A^{p) : si{x,y) — si{y,x)}. 

We will use the sim,plified notation A{p) and As{p) for A'^{p) and A^{p), respectively. 

Definition 2.2. The set of reversible twist maps F of the form (fTp with s G BoiS) C A^{p) will be referred to as 

(4) :F^'Pi~s) = {F:{x,-s{y,x))^{y,s{x,y))\ s € B,{~s) C A^^P)} ■ 

We will also use the notation 

We are now ready to state our main theorems from (Gaidashev et al 2011). The original formulation of Main 
Theorem l2.3l in (Eckmann et al 1984) was in the space yl°-^(1.6). The statement below was proved in (Gaidashev et 
al 2011). 

Theorem 2.3. (Existence and Spectral properties) There exists a polynomial sq : C^ i— > C, such that 

i) The operator TZekw is well-defined, analytic and compact in Bg^^so) C As{p), with 

p = 1.75, go = 5.79833984375 x 10""*. 

ii) There exists a function s* G Br{so) C .4s (p) with 

r = 1.1 X 10"^°, 

such that 

T^EKwls*] = s*. 

iii) The linear operator DTZekw[s*] has two eigenvalues outside of the unit circle: 

8.72021484375 < Si < 8.72216796875, S2 = ^, 

where 

-0.248875313689 < A, < -0.248886108398438. 

iv) The complement of these two eigenvalues in the spectrum is compactly contained in the unit disk. The largest 
eigenvalue in the unit disk is equal to A*, while 

spec{DTZEKw[s*])\{6i,S2,X*}'^ {z £ C : \z\ < 0.1258544921875 = z^}. 

The Main Theorem 12. 31 imply that there exist codimension 2 local stable manifolds yV-jzEKwi^*) ^ .4s(1.75) of the 
operator TZekw- 

Compactness of the operator TZekw in neighborhood of s* implies that there exists a strong "submanifold" 

W^,,„,(s*)cW,z,,.„„(s*), 

of codimension 1 in yViiEKw (**)i such that the contraction rate in W^ , (s*) is bounded from above by the number 
ly from Part iv) from 12.31 

w^'kKwM-T^EKwmp^oiiyn 

for any two s and s in W^ (**)• 

Definition 2.4. The set of reversible twist maps of the form [^ such that s G y^TZEKwi-^*) *- .As(1.75) will be 
denoted Wekw , <md referred to as infinitely renormalizable maps. 

The set of reversible twist maps of the form {!]) such that s G W^ , ,(s*) C .4^(1.75) will be denoted W^^j^yy. 

Definition 2.5. Set, 

WgCs) = WEKwnTl'^^s), 



Q 

1-75 (--A 



Wlis) ^ W%J^^r^T'g'^(s) 
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Naturally, these sets are invariant under renormalization if g is sufficiently small. 

The eigenvector of A, is denoted by ipfJ^"^ = V's* + '0 and hf^^, for DTZekw and DRekw, respectively ( see 
(Gaidashev et al 2011) for details), where tps* is the eigenvector with eigenvalue A, for the operator DTZo induced 
hy Ro[F] = A~^ o F o F o Af. -0 is an eigenvector of DTZq and DTZekw with eigenvalue 1 and 0, respectively. The 
A, direction is constructed by a coordinate change of the form S~^ o F o Se, where 



(5) Se{x,u) = [ X + ex 



„2 



l + 2e2; 



(6) S^\y,v) = [ ^^'~\ vVmry 

The map S^^ o F o Se is generated by: 

s{x, y) = s{x + tx^,y + ey^)(l + 2ey). 

In addition, the following more general fact about the spectrum of the operator 7?.o[s*] holds (see (Gaidashev et 
al 2011)). Denote by 02{'D) the Banach space of maps F : V ^^ C^, analytic on an open simply connected set 
I? C C^, continuous on dV, equipped with a finite max supremum norm || • \\x>'. 

||F||i5 = niax< sup \Fi(x,u)\, sup \F2{x,u)\\ . 

I (2;,u)el' (x,u)£'D \ 

Lemma 2.6. There exists a domain V, such that the operator Rq is differentiable in a neighborhood of F* in 02{'D), 
and spec(-Di?o[^*]) contains eigenvalues /iiA^, * > 0, j > 0, of multiplicity at least 2, and niX^^, j > 0, and X^^^^, 
* ^ 0; of multiplicity at least 1. 

Remark 2.7. Note, that not all of the maps in 02(2?) are generated by generating functions, neither do all of them 
admit the symmetry T o F oT = F^^ , therefore not all of the products of powers of jii, and A* mentioned in Lemma 
\2.6\ are eigenvalues of DRekw[F*] in TF-'J^g''^^{s*). 

The following two theorems describe the dynamics of the infinitely-renormalizable maps on their "stable" Cantor 
sets. 

Theorem 2.8. (Stable Set) 

There exists g > such that any F (E Wq^sq), admits a "stable" Cantor set Cp C T) with the following properties. 

i) For all x G Cp the maximal Lyapunov exponent xi^'i F) exists, is F -invariant, is equal to zero: 

x{x;F)=0, 

and 

1 f \\DF^ix)v\\ \ 
hm — log <^ ^—r \ = 0, 

l-!-±00 \l\ [_ \\V\\ J 

uniformly for all v gM.'^ \ {0} and x G Cp- 
ii) The Hausdorff dimension of Cp satisfies 

diuYH{Cp) < 0.794921875. 

iii) The restriction of the dynamics F\cp is topologically conjugate to the adding machine. 

Theorem 2.9. (Rigidity) Let s* andCp be as in Main Theorems \2.3\ and \8.8\ There exists g> 0, such that for all 
F and F in Wi,{s*), 

where h extends to a neighborhood of Cp as a differentiable transformation, whose derivative Dh is Holder continuous 
with the Holder exponent 

a > 0.0129241943359375. 
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3. H5lder continuity of the derivative 
In this section we will give numerical evidence in support of the following conjecture. 

Conjecture 3.1. For all F,F e Wg„{s*), go ^ 6 x 10"^ 

F\c,lF\c„ 

where h is a dijferentiable transformation, whose derivative Dh is Holder continuous with positive Holder exponent 
larger than: 

a > 0.070347 

We would like to emphasize that it is not the C^+" property that we are after, but rather the specific value of 
a. The fact that the conjugacy extends to a differentiable function whose derivative is Holder has been proved in 
(Gaidashev et al 2011). 

To investigate the Holder exponent of the derivative of the conjugacies in Wgg{so), we conduct an experiment, 
which will be described below. The convergence of Dh^ depends on the spectral properties of the operator DRekw 
in the codimension one submanifold Wg^ , which is transversal to hf^^^ at F* . Note that all our computations are 
done with the entire neighborhood {F S Wg„ ) of F* on the stable manifold. To ensure that we study two different 
functions, we study a F G Wq„ transported along the one dimensional manifold tangent to the eigenvector hf 
at its intersection with W^ , which we denote by I{F). In practice this means that we study conjugacies from an 
arbitrary F e W^^, to an arbitrary F e W^^, transported some distance along I{F). We study such functions using 
the generating function s^, defined by 

U^,y) := s*((l - e)x + e{l - e) V, (1 - e)y + e{l - e)V)(l + 2e(l - e)y)/(l + e) 
= s*{x,y) + exPs* +f^i^ + 0{e^). 

The corresponding scalings are denoted by Ae and /ig, they are computed from the formulae 

(7) s,(A„l) + s,(0,l) = 
and 

(8) Me = Xedis^{Z[s^]{x,y),y)\^^^^^, 
where diZ[s^](x,y) is computed using the midpoint formula: 

Se{x,Z[s^]{x,y)) + Se{y,Z[s^]{x,y)) ==0, 



I.e., 



= di{s^{x,Z[s^]{x,y)) + s^{y,Z[s^]{x,y))) 

= dis^{x,Z[s^]{x,y)) + d2S^{x, Z[se]{x,y))Zi[se]{x,y) 
+ d2S^{y, Z[s^](x, y))Zi[se]{x, y), 



from which 



(9) d^Z[sMK,0)- -^i^^(-^-l) 



follows, and hence 
(10) 



d2S,{Xe,l)+d2s40,l) 

~XAs,{l,0)di.UX,,l) 

d2s4Xe,l)+d2U0,l) ■ 
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To approximate the renormalization operator acting on s^, we note that 

Ti-EKwlse] = TZEKw[s*]+eDTZEKw[s*]ips* + eDTZEKw[s*]'4> + 0{e'^) 
= TlEKwis*] + eDnEKw[s*W + 0(e2) 
= s* + eDno[s*]iJs* + eKi' + 0{e^) 

= s.A. +0(e'). 
Thus, 

= s,(;,.)„+0(6V(l-(A*)')) 

= S,(A.).+0(6'). 

We denote by F^ the map generated by a generating function s^. To siniphfy the notation we denote the corresponding 
presentation functions by 



*0 


:= 


Af. 




n 


:= 


F^ o Af, 






:— 


*^ o • • 





According to (Gaidashev et al 2011), the maps 

converge to a conjugacy of -F*|c. and Ff\cp^, which is, furthermore, extendable to a C^+" map on a neighborhood 
ofC,. 

The foUowing Lemma about the existence of hyperbohc fixed points for maps in a small neighborhood of the 
renormalization fixed point map F* is a restatement of a result from (Gaidashev and Johnson 2009a) in the setting 
of the functional space ^s(1.75). The proof of the Lemma is computer-assisted (see (Gaidashev and Johnson 2009a)). 

Lemma 3.2. Every map F G J^p{so), with g = 6.0 x 10^^^ and p = 1.75, possesses a hyperbolic fixed point p^ G V, 
such that 

1) TTxP^ G (0.577606201171875,0.577629923820496), and tTuP^ — 0, where tTx.u o,re projections on the x and u 
coordinates; 

2) DF{p^) has two negative eigenvalues. 



F 



e I 



G (-2.0576171875, -2.057373046875), 



e^ G (-0.486053466796875, -0.48602294921875), 
corresponding to the following two eigenvectors: 

s^ = [1.0, -(0.77978515625, 0.779815673828125)], and u^ = T(s^). 

This Lemma implies the existence of hyperbolic 2"-th periodic orbits for maps in yVg(so)- Let On{F) denote such 
2"-th periodic orbit of F G Wg{so), specifically: 

2"-l 

0^{F)= U F'«„(p^")), 
where p^" is the fixed point of F„ = _R"[F] G Wg{so). We will also denote 

Consider the 2"-periodic point p^^n = p^„ g -B* „ . 
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We will estimate the Holder exponent of the conjugacy ft,*^ by computing perturbations oi p^n in -B*^. 
To this end, we estimate 

(11) iV„(e,w",t) :== \\Dhl4p^r.) ~ Dhl^p^r. + 6{Lj'')icost,smt))l 
where, for n < 15, 

with 9 = 0.272 turns out to be a good choice. The reason that we multiply with |lpa;"|l is that the size of different 
points on a periodic orbit can vary dramatically, and we want to make the perturbations on one renormalization level 
to be, relatively, of the same order, i.e., we perturb the same significant digit of each point on the periodic orbit. 
We remark that S needs to be sufficiently small to guarantee that the perturbation is inside of -B* „ , and at the same 
time sufficiently large to be computable. In particular we want to be able to compute 7V„(e,aj", t) with n as large 
as possible. Our computations are done with IEEE extended precision (80 bits), for which the 6 chosen above is a 
good compromise, and allows us to compute with n < 15. 

We are interested in bounds on a that are uniform in t and w". Therefore, we define 

(12) M„(e,w") := max iV„(e,w",i) 

te[0,27r] 
.^,^ .^ . l0gM„(£,C.") 

(16) ctnle) := mm — — — . 

^ ' ^^ ^"£{0,1}" log(5(cj") 

The motivation to the definition of a„(e) is that we want to find an a such that 

\\Dhl4x)-Dhl4y)\\<C\\x-yr, 

for all X 7^ y e i?* „ . 

As n ^ oo, (5(a;") -^ uniformly, since diam(i?*„) < C0". Thus, we get the following dichotomy, if 

Q!„(e) -^ a{e) > 0, 

then a exists and is equal to a(e), and if there exists m, such that 

Q!„(e) < 0, for n > m, 

then a does not exist. 

Obviously, we are only able to collect a finite amount of data, i.e., we are only able to compute Nn, for a finite 
discretization of [0, 27r], and only for a finite number of renormalization levels. The first discretization is not a major 
concern, since by discretizing until the results stabilize, we can be almost certain that our results are correct (up to 
some accuracy), this issue will be explained further in the section below. In the results given below we have discretized 
[0, 27r] with 512 points. The computed values of a„(e) agree to 5 digits when the discretization is increased from 256 
to 512 points. It is therefore reasonable to conclude that further discretization of the interval [0, 27r] would have a 
negligible effect on our results. 

The second problem is much greater; ideally we should compute a„(e) for large enough n's so that anie) converges. 
Unfortunately, this is not possible; on one hand due to the enormous complexity of the computation of iV„, which is 
0(2"), and on the other hand due to the finite accuracy of computers. Our choice of 6{uj") = ^^'°^^"||pcj"||: means 
that 6{uj^^) « 3.36 x 10~^^, and since the precision of the 80-bit extended precision is 19 decimals, this is already 
too small to be reliably computable. Note that for some choices of w", p^n., might be macroscopic, regardless of n, 
preventing us from proceeding further. In principle we could use software implemented, higher precision libraries 
to compute for higher n's, but, as we will show below, we would need to compute the first 30 — 50 levels to get 
convergence, and 2^° is a prohibitively large number of combinations. Therefore, we stop at n = 15. 

In order to say something about the convergence of a„(e) from a finite number of renormalization levels, i.e., 
from computing a finite number of /i^„ , we try to fit the results to a graph of a function which approaches a 
constant positive value. It turns out that for any choice of e in the range, where our computation yields results, our 
measurements fit almost perfectly to a function a^ (n) of the form 

2 /^2 I 7,3 /^3 I 7,4 /^4\ 



a,Xn) ~ a.exp (kl/n + kl/n^ + fc^/n^ + k^' 



n 



clearly, a^n) -^ a^. We set 
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This estimation of the measurements by the function ae{n) is illustrated in Figure [Jl where we have restricted to 
the subset 3 < n < 15, the reason being that we would like to avoid transient effects. These transient effects arise 
primarily from the fact that for small n, (5(w") is rather large compared to M„(e,aj"). 

The combined function constructed by gluing together the measurements an(e) and the graph of the function 
ae(n) at n = 15 is illustrated in Figure [TJ From Figure [T] we see that the function does not flatten out completely (up 
to the resolution of the image) until n — 30; and as already mentioned, computing /i^3o, is completely unrealistic. 
The computation at our final rcnormalization level, n — 15, already takes about 4 hours for one choice of e. 



Figure 1. (a)Tracking the measurements: the blue curve is the measured value of Q;„(e) as a 
function of the renormalization level. The green curve is the estimate of the form a^{n) = 
fle exp (k/n + kl/n? + k^/n? + k^/rr^. The results in this figure are from the computations with 
e = 1.593584796420859 x 10^*. (b) Extrapolation of the measurements: to extrapolate the mea- 
surements, the measured data, a„(e), is glued together with the graph of the computed estimate, 
a^{n). The blue curve is the measured value of a„(e) as a function of the renormalization level, the 
red star is the last measurement. After the red star, the green curve is the graph of the function 
a^{n) — Qe exp (fc/n + kl/n"^ + k^/n^ + k^/nr^. The results in this figure are from the computations 
with e = 1.593584796420859 x 10"^ 

We implement an algorithm based on the discussion above in a C-I--I- program, (ProgsB 2011), and compute the 
estimated rigidity of /i^ for the largest possible range of e. This range turns out to be [10~^, 10~^]. The result of the 
computations together with the relative least squares error of the estimation is given in Table [T] 



e 


a{e) 


13 V ^"=3 


a„(e) 

On 


-a^(n) 12 


10"'* 


3.37 X 10-^ 






0.40% 


10-5 


1.13 X 10-2 






0.25% 


IQ-^ 


2.21 X 10-2 






0.19% 


10-7 


3.59 X 10-2 






0.13% 


1.593584796420859 x 10"^ 


7.03 X 10-2 






0.50% 



Table 1. The estimated values of a(e) for various values of e, together with the least squares error 
of the estimating function a^in). 



To discuss how the rigidity of the conjugacies between F^ and F* depends on the distance of F^ from the fix point; 
we need to relate the size of e with ||se ~ s*||- This correspondence can be estimated using the fact that: 
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^iiJMiA^^m^ 



Figure 2. Distribution of angles in Oi2{F*) 



Thus, 

By using the estimate 



.*||=e||V^--^||+0(e^) 



llV-f.^^ll > 37.6509616148184234, 

together with the estimate given in Table [T] we conclude that go corresponds to an e < 1.593584796420859 x 10^* 
and that 



F' 



This conclusion implies Conjecture 13.1 



\cpe , F*\c,, ioT F'- eIcWe„{sQ), with a = 0.070347113466231. 



4. Twisting of line fields 



Consider the stable and unstable invariant direction fields on the 2"-th periodic orbit On{F). At every point p^, 
00 e {0, 1}" of On{F), these directions are given by 

,F/F„\„F„ 



(14) 
(15) 



u 



= i?<(p^")u^". 

The angles between these vectors and the positive real line will be denoted by a^ and /?^, respectively. 
We have plotted the distribution of the angles a^ for lo G {0, 1}^^ for vectors (fT4| in Fig. [2l 
The figure demonstrates that the stable (and, via reversibility, the unstable) direction undergoes a sort of twisting: 
the directions do not concentrate, but are distributed with a certain measure over the whole circle. 
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One can heuristically explain the structure of this distribution. First, the appearance of two large peaks at angles 
and TT can be explained by the fact that all vectors s^ with in the beginning of their coding u tend to be aligned 
horizontally: indeed, for every w € {0, l}""'"', 

fe 



tana, 



Qfc 



const 



tana 



Fk-i 



The smaller peaks (such as the two peaks between angles 1 and 2) correspond to the stable directions at points p[^ 
that are direct and inverse images under dynamics F of points p^ whose angles are in one of two central peaks. 

Let a V be a unit vector. We identify the set of unit vectors in M^ with T^, and write w £ T^. Consider the 
normalized derivative cocycle ^^ : Ci? x T^ !—> T^: 

DF\x)v 

Ai{x)v= ^ ' 



\DF^{x)v\\ 



Fig. [2] lends evidence to the following conjecture: 



Conjecture 4.1. There exists an absolutely continuous invariant measure v onT^ , which is ergodic for the cocycle 

Ai{x), X e Cf-' 

n-l 



(16) 



lim i^/(A(a;»^^/ 



f{y)dv{y), 



for V a. a. points v Q T^ and every v-measurahle and integrable / : T"'^ i— > M. 

To support the conjecture we have computed the left and right hand sides of the equation (J16p for several monomials 
in cos and sin. The approximation of the density of the measure has been computed by dividing the circles into a 
collection of iV subintervals [9i, 9i+i] of equal length, and counting the relative number K" of angles a^ , to € {0, 1}^^, 
in On that lie in those intervals. 

We have collected the data in Tables [Sl—TSl The number reported in the tables is |2(L — R)/{L + R)\, where 

Af-l N-l 



^ = ]^E/(-^fc«0,0))^;), R^J2f 



fc=0 



i=0 



e, + e. 



4+1 



K?. 



Here, we have used that the point (0, 0) is in the stable Cantor set (see (Gaidashev and Johnson 2009b) for details). 
The quantity L has been computed for two different w's. 



Or, 



71= 12 

n = 14 

71= 16 



N = 1000 
2.5767 X 10-^ 
2.0387 X 10-3 
2.5675 X 10-4 



N = 5000 
3.0478 X 10-3 
2.0535 X 10-3 
9.8508 X 10-'^ 



N 



15000 

3.0450 X 10-3 
2.0552 X 10-3 
1.9443 X lO-'' 



Table 2. The relative difference of the left and right hand sides in ([T6l) for the function sin^, 
M = 20000, u= (1,0). 



Or 



n = U 

71= 16 



N = 1000 
4.3071 X 10-3 
2.9367 X 10-3 
2.4732 X 10-4 



N = 5000 



4.2824 X 10-3 
2.9642 X 10-3 
2.1161 X 10-4 



N = 15000 



4.2809 X 10-3 
2.9731 X 10-3 
2.0896 X 10-4 



Table 3. The relative difference of the left and right hand sides in ([T 
M = 20000, u= (1,0). 



for the function cos 
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On 



n = 12 
n = 14 

71= 16 



N = 1000 
2.6290 X 10-^ 
1.6695 X 10-3 
1.1324 X 10-4 



iV = 5000 
3.1001 X 10-3 
1.6836 X 10-3 
4.6830 X 10-* 



N = 15000 
3.0973 X 10-3 
1.6854 X 10-3 
5.6424 X 10-"* 



Table 4. The relative difference of the left and right hand sides in ([T 
M = 20000, V = (0, 1). 



for the function sin 



On 



n = 12 
n = 14 
n = 16 



N = 1000 
4.4239 X 10-3 
2.0916 X 10-3 
1.0927 X 10-3 



N = 5000 



4.3992 X 10-3 
2.1189 X 10-3 
1.0568 X 10-3 



N ^ 15000 



4.3978 X 10-3 
2.1278 X 10-3 
1.0542 X 10-3 



Table 5. The relative difference of the left and right hand sides in p6 
Af = 20000, w= (0,1). 



for the function cos 



5. Reality of the renormalization spectrum 

It has been a long standing conjecture, supported by numerical evidence, that the renormalization spectrum of 
the Fcigcnbaum period-doubling operator for unimodal maps is real. We observe the same phenomenon with the 
period-doubling renormalization spectrum for area-preserving maps 

We recall that our Theorem 12.31 has been proved for the operator TZekw acting on the space ^s(l-75). 

We will choose a new basis {i^tj} in As{p)- Given s € As{p) we write its Taylor expansion in the form 






where ipij € As{p)- 






-1, i>o, 



x'+^y^, i 

x'+'y'+'-^x^+'y\ *>-l, j>i, 
J + 1 



i,. 



V'»,j 



IJ 



■, i > -1, j>max{0,i}. 



and the index set / of these basis vectors is defined as 

I = {(i,j) e Z^ : i>-l, j > niax{0,i}}. 
We will denote the cardinality of the set 

lN = {{i,j)eI:i+j <N} 
as D{N). Assign a single index to vectors ipij, {hi) G Ini as follows: 

fc(-l,0) = l, fc(-l,l) = 2, ..., fc(-l,A^) = A^ + 1, fc(0,0) = A^ + 2, 



fc(0, l) = 7V + 3, 



N -1 



,N-l- 



N - 1 



= D{N). 



This correspondence («, j) i— >■ fc is one-to-one, we will, therefore, also use the notation {i{k),j{k)). 

For any s G As{p), we define the following projections on the subspaces of the linear subspace Ej)(^^^ spanned by 

r ; l-D(W) 

m<D{N) 
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We have found an approximation of a matrix representation of the finite-dimensional hnear operator 



Hf 



DTI 



EKW 



:*inso( 



as 



Dn,m = ll-mDTZEKw[s*]ll^n, 

and have computed the eigenvalues of the matrix. Table[8]gives the first 100 eigenvalues of the o-pevsAoi He o{n)^'^ekw[s*]^e ui 
N = 20. 

We would like to note that the complexity of some of the higher-order eigenvalues is likely a numerical artifact. 
To demonstrate this, we list eigenvalues 33 and 34 of the operator IIej^^j^-iDTZekw[s*]Hed{n) fo^ several values of N 
in Table [6l and eigenvalues 62 and 63 in Table [7l One can see that as one increases the dimension of the projection, 
the pair of complex conjugate eigenvalues turns real after some iV, and stays that way. 



N 


33 


34 


10 


7.710128558752339 x 10"' - 1.061214359986916 x 10" 'i 


7.710128558752339 x lO"' + 1.061214359986916 x lO-'i 


15 


9.116707201322168 x 10"' +0.0i 


8.866211124768190 x 10"' + O.Oi 


20 


9.116282429256495 x 10"' + O.Oi 


8.873952087973409 x 10"' + O.Oi 



Table 6. A comparison of the eigenvalues 33 and 34 of the operator IIejj,,^,DTZekw[s*]^Eoin) ^^^ 
several values of N. 



N 


62 


63 


15 


-6.705013355269772 x 10""" - 1.625832021663626 x 10"'"i 


-6.705013355269772 x 10""" -1- 1.625832021663626 x 10-""i 


20 


-8.704206290836820 x 10""' + O.Oi 


-8.227051108116014 x 10""" + O.Oi 


25 


-8.704155784717757 x 10""" + O.Oi 


-8.247546278128800 x 10""" + O.Oi 



Table 7. A comparison of the eigenvalues 62 and 63 of the operator IIejj,j^^D'R'EKw[s*]TIeoin) ^'^^ 
several values of N. 
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1 


8.72109720060341027 + O.Oi 


iV 


51 


-1.36773104970236911 x 10"** + O.Oj 






Y 


2 


-4.01807670479890925 + O.Oi 


Y 


52 


-1.29617185901792936 x 10"* + O.Oj 






Y 


3 


-2.48875288718523027 x 10^^ + O.Oi 


Y 


53 


-1.33177097864396282 x 10"** + 0.0? 






Y 


4 


-1.16629420927308277 x 10^^ + O.Oi 


N 


54 


5.83632841811164741 x 10"'' + O.Oj 






N 


5 


7.29842918134375602 x 10^^ + O.Oi 


N 


55 


3.49734698280496765 x 10"'' + 0.0? 






Y 


6 


6.19389093347281188 x IQ-^ + o.Oi 


Y 


56 


3.40553232521175106 x lO"'-* + O.Oj 






N 


7 


-1.54150639435907658 x 10^^ + O.Oi 


Y 


57 


3.30875856445169904 x lO"'^ + 0.0? 






N 


8 


-1.50053025013312190 x 10"^ + O.Oi 


Y 


58 


3.23415716696474633 x lO"'* + O.Oj 






N 


9 


-5.24316259408011873 x 10"^ + O.Oi 


N 


59 


3.20031850876952441 x 10"'' + 0.0? 






N 


10 


3.83642848957628521 x 10"^ + O.Oi 


Y 


60 


-1.21497298394740393 x 10"'' + O.Oi 




N 


11 


3.73444899232657950 x 10"^ + O.Oi 


Y 


61 


-8.47423932266715023 x 10"''' + O.Oi 




Y 


12 


1.67255839935358406 x 10"^ + O.Oi 


N 


62 


-8.70420629083682033 x lO""' + O.Oi 




Y 


13 


1.95630639839669629 x 10"^ + O.Oi 


N 


63 


-8.22705110811601487 x 10""' + O.Oi 




N 


14 


-9.54792247990135266 x IQ-* + O.Oi 


Y 


64 


-8.07179147279961451 x 10""' + O.Oi 




N 


15 


-9.29412071169940800 x IQ-* + O.Oi 


Y 


65 


-7.79143693124703818 x 10""' + O.Oi 




N 


16 


-7.74683770828318898 x 10"* + O.Oi 


N 


66 


4.63752609505868076 x 10"'" + O.Oi 




N 


17 


2.37624196383790100 x 10""* + O.Oi 


Y 


67 


-3.22686335036190861 x 10"'" + O.Oi 




N 


18 


2.31307697553682750 x lO"* + O.Oi 


Y 


68 


2.77581537948968217 x 10""' + O.Oi 




N 


19 


2.05791964161655016 x 10""* + O.Oi 


N 


69 


2.16643085621470110 x 10""' + O.Oi 




Y 


20 


4.61021338856759535 x lO"'"'' + O.Oi 


N 


70 


2.09856849656907084 x 10""' + O.Oi 




N 


21 


-5.91387904855386636 x 10"^ + O.Oi 


Y 


71 


2.08109796932933800 x 10"'" + O.Oi 




N 


22 


-5.75667700266775337 x 10"^ + O.Oi 


Y 


72 


1.96066223947378410 x 10"'" - 4.43072486498362171 


X 10"'2i 


N 


23 


-5.60365367884625166 x 10"^ + O.Oi 


Y 


73 


1.96066223947378410 x 10"'" + 4.43072486498362171 


X 10"'2i 


N 


24 


-3.87086314650991711 x 10"^ + O.Oi 


N 


74 


-4.00724308724832731 x 10"" - 1.29484585578379769 


X 10""i 


N 


25 


1.47181835587117418 x 10"^ + O.Oi 


Y 


75 


-4.00724308724832731 x 10"" + 1.29484585578379769 


x 10""i 


N 


26 


1.43269464891488572 x 10"^ + O.Oi 


Y 


76 


-5.39409517510675544 x 10"" + O.Oi 




Y 


27 


1.39461092667464458 x 10"^ + O.Oi 


Y 


77 


-5.22615941114597117 x 10"" - 6.44631270712982928 


x 10"'^i 


N 


28 


3.37032628919152810 x 10"^ + O.Oi 


N 


78 


-5.22615941114597117 x 10"" +6.44631270712982928 


x 10"'-*i 


N 


29 


-3.66299218723578256 x 10"® + O.Oi 


Y 


79 


-4.83301369903105249 x 10"" - 1.61011502564644013 


x 10-'"i 


N 


30 


-3.56562288094160436 x 10"® + O.Oi 


Y 


80 


-4.83301369903105249 x 10"" + 1.61011502564644013 


X lO-'^i 


N 


31 


-3.47084200041532543 x 10"® + O.Oi 


Y 


81 


-1.21652373238263412 x 10"" + O.Oi 




N 


32 


-2.68212102338467744 x 10"® + O.Oi 


N 


82 


1.34403330160617114 x 10"" + O.Oi 




Y 


33 


9.11628242925649478 x 10"^ + O.Oi 


Y 


83 


1.29886814444646485 x 10"" - 4.49484074640346900 


x 10"'^i 


N 


34 


8.87395208797340851 x 10"'' + O.Oi 


Y 


84 


1.29886814444646485 x 10"" + 4.49484074640346900 


X 10"'3i 


N 


35 


8.63807486025034232 x 10"^ + O.Oi 


Y 


85 


1.30874465804588329 x 10"" + O.Oi 




N 


36 


8.43743910455452052 x 10"^ + O.Oi 


N 


86 


1.22553839952396333 x 10"" + O.Oi 




N 


37 


-6.13340779558256008 x 10"^ + O.Oi 


N 


87 


1.14128486730559826 x 10"" - 1.22675413039504159 


X 10"'^i 


N 


38 


3.64333304618918514 x 10"^ + O.Oi 


N 


88 


1.14128486730559826 x 10"" + 1.22675413039504159 


X 10"'^i 


N 


39 


-2.26881741829139139 x 10"^ + O.Oi 


Y 


89 


-4.64857574616255465 x 10"'^ - 4.15935498398141420 


X lO-'^i 


N 


40 


-2.20850766682014296 x 10"^ + O.Oi 


Y 


90 


-4.64857574616255465 x 10"'^ + 4.15935498398141420 


X 10"'2i 


N 


41 


-2.14980452161630267 x 10"^ + O.Oi 


Y 


91 


-3.50139347550126033 x 10"'^ + O.Oi 




N 


42 


-2.09265346840381688 x 10"^ + O.Oi 


Y 


92 


-3.12844287139405750 x 10"'^ - 2.35599032104882566 


X 10"'3i 


Y 


43 


5.64652471456064890 x 10"** + O.Oi 


Y 


93 


-3.12844287139405750 x 10"'^ + 2.35599032104882566 


X 10"'^i 


Y 


44 


5.49645085465527861 x 10"** + O.Oi 


Y 


94 


-3.18763936883361579 x 10"'^ + O.Oi 




N 


45 


5.35030991885504657 x 10"** + O.Oi 


Y 


95 


-2.60876535503609472 x 10"'^ + O.Oi 




N 


46 


5.20806572338405750 x 10"** + O.Oi 


Y 


96 


1.00716118446860206 x 10"'^ + O.Oi 




N 


47 


4.06666476717650202 x 10"** + O.Oi 


N 


97 


-8.13328302601060092 x 10"'^ + O.Oi 




N 


48 


-2.89661341440091808 x 10"** + O.Oi 


N 


98 


8.04488862233887989 x 10"'^ - 1.00262940606313244 


X 10"'^i 


N 


49 


-1.44725130238278229 x 10"** + O.Oi 


N 


99 


8.04488862233887989 x 10"'^ + 1.00262940606313244 


X 10"'^i 


N 


50 


-1.40530550882139069 x 10"** + O.Oi 


Y 


100 


5.08237476178472880 x 10"'^ + O.Oi 




Y 



.DUf 



iHf 



Table 8. Some of the eigenvalues of IIei^^i^^-L>Kekw[s \i^i-Eo(N)^ 

those eigenvalues that correspond to the products of powers of A* and /i* . 



N = 20. The last column hsts 
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